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We investigate theoretically BEC-BCS crossover physics in the presence of a Rashba spin-orbit 
coupling in a system of two-component Fermi gas with and without a Zeeman field that breaks 
the population balance between the two components. A new bound state (Rashba pair) emerges 
because of the spin-orbit interaction. We study the properties of Rashba pairs using a standard pair 
fluctuation theory. At zero temperature, the Rashba pairs condense into a macroscopic mixed spin 
state. We discuss in detail the experimental signatures for observing the condensation of Rashba 
pairs by calculating various physical observables which characterize the properties of the system and 
can be measured in experiment. 

PACS numbers: 05.30.Fk, 03.75.Hh, 03.75.Ss, 67.85.-d 



I. INTRODUCTION 

Since its recent realization in cold atomic systems [1|- 
artificial gauge field has received tremendous attention. 
The concept of gauge field is ubiquitous, a classical ex- 
ample of which is electromagnetism. The achievement of 
the above mentioned experiments allows us to simulate 
charged particles moving in electromagnetic fields using 
neutral atoms. The more recent realization of a non- 
Abelian gauge field in a system of 87 Rb condensate [3] 
provides us a system of spinor quantum gas whose inter- 
nal (pseudo-)spin degrees of freedom and external spatial 
degrees of freedom are intimately coupled. Novel quan- 
tum states will emerge in such spin-orbit coupled sys- 
tems [fj(. Although experiments on artificial gauge field 
have so far only been carried out in bosonic systems, we 
have no reason to doubt that they will soon be extended 
to fermionic systems. Theoretically, there have been a 
number of papers focusing on the interesting properties 
of spin-orbit coupled Fermi gases in both three dimension 
(3D) [Ml and 2D [20H23j. most of the them appeared 
just over the past few months. 

The salient features of spin-orbit coupled fermions in- 
clude: enhanced pairing field [H H, [Hj], mixed spin pair- 
ing [24| , non-trivial topological orde r fl2l . [25j , and pos- 
sible existence of Majorana fermion [26|, etc. The pur- 
ose of the present paper is to extend our earlier work 
by providing a detailed description of the theoretical 
techniques and by including the effect of a Zeeman field 
which not only breaks the population balance, but also 
may induce topological phase transitions in the system. 
We start from a discussion of the two-body problem, fol- 
lowed by a detailed study of the many-body system. We 
present our calculations of various important physical ob- 
servables such as the single-particle spectrum, density of 
states, spin structure factors, etc., which may be used to 
characterize the system experimentally. 

The paper is organized as follows. The model Hamil- 
tonian and the general technique are presented in Sec. ITT] 



The results for the two-body and the many-body prob- 
lem are presented in Sec. lIIII and llVl respectively. Finally 
we provide an outlook and conclusion in Sec. [V] 



II. PHYSICAL MODEL AND GENERAL 
TECHNIQUE 

In this section, we first present the model Hamiltonian 
of interest and then give a detailed description of the 
functional path integral formalism employed in deriving 
the relevant equations. We choose this formalism as it 
allows us to present a unified treatment for both the two- 
body and the many-body physics. 



A. Model Hamiltonian 

Here we consider the BEC-BCS crossover theory in the 
presence of a Rashba spin-orbit (SO) coupling X(k y a x — 
k x a y ) in the x-y plane, together with an external Zee- 
man field ha z . The Zeeman field acts to break the pop- 
ulation balance between the two spin components of the 
fermions. The second-quantized Hamiltonian for a uni- 
form system reads, 



H 



dr- 



£k + ha z + \{k y a x - k x a y ) 



1> 



+ U^ J +{Y)i,+ {v)i, i {v)^{v)Y 



(1) 



where £k = h 2 k / (2m) — \i with /z being the chemi- 
cal potential, and tp (r) = (r) , ip± (r)] T , \jj a (r) is the 
fermionic annihilation operator for spin-cr atom. Here h 
is the strength of the Zeeman field and A is the Rashba 
SO coupling constant. Without loss of generality, we 
take both h and A to be non-negative. The last term in 
Eq. (JTJ represents the two-body contact s-wave interac- 
tion between un-like spins. 
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The SO coupling term has already some interesting 
effect on the single-particle physics. The single-particle 
spectrum (i.e., the eigenenergy of the dressed states) can 
be straightforwardly obtained as 



the ground state becomes non-degenerate and occurs at 
momentum k = 0. 



e k ,± ^^±Jh 2 + X 2 k 2 ± 



(2) 



B. Functional Path Integral Formalism 



where fc^ = J k 2 + k 2 is the magnitude of the transverse 
momentum. The lowest single-particle state occurs at 



and 



^/m 2 X 2 /h A - h 2 /X 2 , h<m\ 2 /h 2 



0. 



otherwise 



(3) 



with the corresponding lowest single-particle energy as 
(taking fi = 0): 

_ f -m\ 2 /{2h 2 ) - h 2 h 2 /{2mX 2 ) , h < mX 2 /h 2 



-h, 



otherwise 



(4) 

Hence for h smaller than a threshold value mX 2 /h 2 , 
the single-particle ground state is infinitely degenerate 
and occurs along a ring in momentum space centered at 
k = and lies in the transverse plane. The radius of 
the 'Rashba ring' decreases as h is increased and van- 
ishes when h exceeds the threshold value, in which case 



We employ the functional path integral method [27, 28] 
to study the problem and consider the partition function, 



Z = \V\i> (r, r) , 4, (r, r)] exp {-S [ff> (r, r) , ^ (r, r)] } , 

(5) 



where the action 



Sty,® = / dT 

Jo 



dr ^2 i>° ( r > T ) 9 T ^a (r, r) + U (if), ip) 



is written as an integral over imaginary time r. Here 
[3 = l/(kBT) is the inverse temperature and % is 
obtained by replacing the field operators %p + and ip with 
the grassmann variables tp and ip, respectively. We can 
use the Hubbard-Stratonovich transformation to trans- 
form the quartic interaction term into the quadratic form 
as: 



J 



e -u f dxdri,^^^ = y x> [A, A] exp J J dr J dr 



|A(r,r)r 
U 



(6) 



from which the pairing field A (r, r) is defined. 

Let us now formally introduce the 4-dimensional Nambu spinor $ (r,r) = [tp-f, tp^,ip-f, ip^] and rewrite the action 
as, 



Z = / £>[$, $; A, A] exp <^ - dr dr dr' / dr' 



-5{r - r')S( T - t') 



where V is the quantization volume and the single-particle Green function is given by, 



(7) 



g- l = 



-9 T - Ck - ha z - X{k y d x - k x a y 
—iAa y 



iAa y 

-d T + + ha z - X(k y a x + k x a y ) 



S(r - t>)6(t - t>) , (8) 



where the Pauli matrix &i (i = 0,x,y,z) describes the 
spin degrees of freedom. Here the momentum k a (a = 
x, y, z) should be regarded as the operators in real space. 

Integrating out the original fermionic fields, we may 
rewrite the partition function as 



Z= \ P[A,A] exp{-5 cff [A, A]} , 



where the effective action is given by 




where the trace is over all the spin, spatial, and temporal 
degrees of freedom. To proceed, we restrict to the gaus- 
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sian fluctuation and expand A (r, r) = Ao + <5A(r, r). 
The effective action is then decomposed accordingly as 
S e s = Sq + AS, where the saddle-point action is 

s °=iV*(-f)->°[- s °"] + £i>. 

(9) 

where Q 1 has the same form as Q 1 in Eq. ((5]) with A 
replaced by Ao, and the fluctuating action takes the form 



with 



E = 



iSAa y 
-i5Aa y 



being the self energy. 



C. Vertex Function 

The low-energy effective two-body interaction is char- 
acterized by the vertex function, which we derive in this 
section. We shall consider the normal state where the 
pairing field vanishes, i.e., Ao = 0. In this case, the in- 
verse Green function Qq 1 has a diagonal form and can 
be easily inverted to give : 



Go(k) 



g+(k) 
5_(fc) 



(10) where k = (k, iuj m ) and 



(11) 



9+{k) = 

g-(k) = 



i 



_ iUm — 6c + ha z + \{k y & x - k x <7 y ) 

iu m - £k - ha z - \(k y a x - k x a y ) (i Wm _ £ k ) 2 - [ft2 + \2 (^2 + £2)] 

1 _ iu m + 6c - ha z + X(ky(7 x + k x a y ) 

+ 6c + ha z - \(k y a x + k x u y ) ~ (iuj m + Z k ) 2 - [/i 2 + A 2 + k{)] 

I 



tui. 



(12) 
(13) 



After some algebra, we may obtain the fluctuating part where the inverse vertex function is given by 
of the action as 



^S = k B Ty Yl [-r-'(9)]*A( 9 )JA(9), (14) 



1/2 



1/2 



(iu) m - £k,+) {iv n - lUJ m - e q -k,+) (i^m - e k~) fan ~ iu m - £q-k,-) 



where ek,± are the single-particle spectrum in Eq. @ and 



A - 



y/W+X^klJh 2 + A 2 (q - kY ± + h 2 + X 2 k x (q x - k x ) + X 2 k y {q y - k y ) 



r- 1 (?) 



47rft 2 a s 21/ 



(iw m - £k, + ) (icj m - £k,-) [}Vn - U*l m - £q-k, + ) fan - iuJ m ~ e q-k,-) 

The summation over iui m in Eq. (|15|) can be done explicitly, after which we find that, 

/ (e q /2+ k ,+) + / (e q / 2 -k,+) - 1 . / (e q /2+k,-) + / (e q / 2 -k ,-) - 1 1 



--E 



W n — e q / 2 +k,+ ~~ e q/2-k,+ 

h 2 + \ 2 (q 2 jA-kl) 



h 2 + A 2 (q/2 + k) \Jh? + A 2 (q/2 - k)* 



w n — e q/2+k,- — e q/2-k, 

C res (q, ii/„;k), 



where fix) = \ j(eP x + 1) is the Fermi distribution function and 

s | [/ ( £ g/2+k,+) + / (e q /2-k,+) - 1] | [/ (fq/2+k,-) + / (e q /2-k,-) - 1] 

W n — £ q /2+k,+ — e q/2-k.+ W n ~ e q/2+k,- — e q/2-k,- 

[/ (eq/2 + k, + ) + / ( £ q/2-k,-) - 1] [/ (Cq/2+k,-) + / ( £ q/2-k, + ) ~ l] 



W n — £ q /2+k : + — e q/2-k- 



W n — £q/2+k,- — e q/2-k,+ 



(15) 



(16) 



(17) 



(18) 
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In writing the above equations, we have replaced the bare 
interaction strength Uq in favor of the s-wave scattering 
length a 8 using 

1 to 1 x - 1 

k 

with e k = h 2 k 2 /(2m). 



III. RESULTS ON TWO-BODY PROBLEM 



Let us first consider the two-body problem. The cor- 
responding two-body inverse vertex function can be ob- 
tained from Eq. (JTTJ) by discarding the Fermi distribution 
function and by setting chemical potential fi = 0. This 
leads to 



47rft 2 a, 2V 



1 



1 



E 



W n — E q /2+k,+ — e q/2-k,+ w n — e q/2+k,- — e q/2-k, 

h 2 + X 2 (ql/A-kl) 



+ 



£k 



h 2 + X 2 (q/2 + k) 2 ±x h 2 + X 2 (q/2 - k)* 



C£(q,u/ n ;k), 



(19) 



where 



■26 

res 



1 



1 



— e q/2+k,+ — e q/2-k, 
1 



W n — e q /2+k,- — e q/2-k,- 
1 



iv n - e. 



q/2+k,+ - e q/2-k,- 



e q/2+k,- — e q/2-k, 



(20) 



One important question concerning the two-body system is whether there exist bound states. The zero- momentum 
bound state energy Eb can be determined from the vertex function using the following relation (iv n — > lu + i0 + ): 



Re[r 2b 1 (q = 0;w = E B ) = 0] 
from which we may derive the equation for the bound state energy as 



47rft 2 a, 2V 



•v ^ 



1 



Eb — 2ek.+ Eb — 2ek,- tk 



1 E' 



Ah 2 



V ^ (E B - 2e k ) (E B - 2e k ,+) (E B - 2e k ,_) 



(21) 



0. (22) 



A bound state exists if its energy satisfies 
Eb < 2e m in , 

where e m in is the lowest single-particle energy defined in 
Eq. Q. 

We solve Eq. (f2"2"j) numerically to find Eb and the re- 
sults are shown in the left panel of Fig. [TJ where we plot 
Eb as a function of the contact interaction strength for 
three different values of the Zeeman field (h = 0.5, 1.0 
and 2.0 toA 2 /h 2 ). For h < mX 2 /h 2 , i.e., when the Rashba 
ring exists as the lowest single-particle state [see Eq. <j3j) ] , 
we always find one bound state solution regardless of the 
sign of a s . It is well known that in the absence of the 
SO coupling, two-body bound state does not exist on the 
BCS side (i.e., a s < 0). The existence of the Rashba ring 
induced by the SO coupling enhances the density of states 
near the single-particle ground state and favors the for- 
mation of bound state By contrast, for h > mX 2 /h 2 , 
the Rashba ring collapses to a point and two-body bound 
state only occurs on the BEC side. Furthermore, the 



larger the h is, the stronger attractive interaction (i.e., 
larger aj 1 ) is required to have a bound state. For ex- 
ample, at h = 1.0 toA 2 /H 2 , the bound state exists for 
aj 1 > 0; at h = 2.0 mX 2 /ti 2 , the bound state exist only 
for a' 1 > 1.35 mX/h 2 (see Fig. [J). 

If there is a tightly bound state, the inverse vertex 
function is simply the Green function of the bound pair. 
A bound state can therefore be examined by calculating 
the phase shift [H, HI 

*(q,a;) = -Im{ln[-r^ b 1 (q ! ii/„^-a; + iO + )]} . (23) 

In the right panel of Fig. [TJ we display the phase shift 
at q = 0. When a bound state occurs, the phase shift 
will have a discontinuous jump of tt at the correspond- 
ing energy as can be seen in the figure. The phase shift 
calculation via Eq. (|2"31 and the bound state energy cal- 
culation via Eq. (|2"2"|) thus corroborate with each other. 

An important quantity that characterizes the proper- 
ties of the bound state is its effective mass. At small mo- 
mentum |q|, we may assume that the bound state has a 
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FIG. 1: (color online). Left panel: Bound state energies Eb as 
functions of scattering length for different Zeeman field h. The 
horizontal dashed lines represent the threshold energy 2e m i n . 
Right panel: Corresponding two-body phase shifts <5(q = 0, uS) 
for different Zeeman field h. Eb and h axe in units of mA 2 /S 2 , 
a s is in units of h 2 /(m\). 




o 
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FIG. 2: (color online). Effective mass 7 = M±/(2m) as func- 
tions of scattering length for different Zeeman field h. h is in 
units of mX 2 /h 2 , and a s is in units of h 2 /(m\). For h > 1, 
the two-body bound state only exists for a s > 0. 



well-defined dispersion, eq = ft 2 g^/(2Mj_) + h 2 q 2 J{4m), 
where Mj_ is the effective mass in the transverse plane. 
Due to the nature of the Rashba SO coupling, the ef- 
fective mass of the bound state along the z-axis is sim- 
ple twice of the atomic mass and is not affected by 
the spin-orbit term. For a given q, we determine Cq 
from the equation: Re [T^ -1 (q; lu = Eb + £q ) = 0] .By 
Taylor expanding the two-body inverse vertex function 
around |q| = 0, and after some tedious but straightfor- 
ward calculation, we obtain: 



1 

7 



2m 
Ml 



4mA 2 Y 



(24) 



where 



X = 



Y 




X 2 ki 



h 2 



\ 2 k\ 



1 



1 

- 2e k ,+) 
1 



1 

{E B - 2e k ,- 



h 2 \ 2 k 2 ± 



{E B - 2e k ,+) {E B - 2e k ,_) {E B 
{E B - 2e k ) 



2e k ) {h 2 + X 2 kl) 2 {E B ~ 2e k )- 



2e k ,+) 2 (E B ~ 2e k ,_) 



Figure [2] displays the effective mass M± as functions of the scattering length a s for several values of the Zeeman 
field strength {h). M± monotonically decreases as l/a s increases. In the BEC limit where a s — > + , M± — > 2m 
independent of the value of h. 



IV. RESULTS ON MANY-BODY PROBLEM 



We now turn to the discussion of the many-body properties. In this work, we only consider the mean-field properties 
of the system, while the effect of fluctuations will be studied in the future. In the mean-field level, the order parameter 
A is a constant and the corresponding momentum space single-particle Green function takes the form (fc = k, iuj m ), 

iu m - £ k - ha z ~ \{k y a x - k x a y ) iA a a y 

-iA Q a y iu m + Ck + ha z - \{k y a x + k x a y ) 



Go 1 (*) = 



6 



Plugging this into Eq. ©, we may obtain the thermodynamic potential as 



,-E ka /k B T 



(26) 



k,a=± 



where E-^± — y + A 2 , + ft, 2 + X 2 k\ ± 2y (ft 2 + A 2 fc^)^ + ft 2 A 2 , is the quasi- particle dispersion. The chemical po- 
tential and order parameter should be determined by, 



= 



N = - 



dfl 



dA ' d[i 
from which we derive the gap and the number equations as follows: 



- - 1 E 



2/(.Eka) - 1 



k.a— ± 



1 + a- 



v/(ft 2 + >?MM + h2A 'o 



N 



V V 



a=± 



K 1 + \ 2 k 



2J.2 



(27) 
(28) 



In the following, we will show various quantities of 
physical interest, obtained from solving Eqs. (|2"7|> and 
([28]) self-consistently. We focus on the zero-temperature 
case, although Eqs. (|27p and (|2"5|) are valid for finite tem- 
peratures. 

We show in Fig. [3fa) and (b) the chemical potential 
\x and the pairing gap A , respectively, as functions of 
the scattering length for different values of the Zeeman 
field. The Zeeman field tends to suppress the pairing 
gap. In Fig. [3jc), we plot the population of the spin- 
up component. For zero Zeeman field, we always have 
equal population in both spin components. For h > 0, 
spin-up component has less population. In all cases, the 
effect of the Zeeman field reduces as the BEC limit (i.e., 
l/a s — > +oo) is approached. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



V(aie) l/(a«fc) l/(a.,fc) 

FIG. 3: (color online). Chemical potential /i (a), pairing gap 
Ao (b), and population of spin-up component rif (c) as func- 
tions of scattering length a s for different values of the Zee- 
man field h at XUf/Ef = 2. Here fcjr = (S^n) 1 / 3 and 
E F = h 2 k 2 F /{2m) are the Fermi momentum and Fermi en- 
ergy, respectively. 

The quasi-particle dispersion E^± is plotted in Fig. [4] 
The spectrum is sensitive to the polar angle 8 of the 
momentum vector k. For k along the z-axis (i.e., 9 = 
and =0), .Ek- may become zero at certain values of k 
when h is sufficiently large, signaling a gapless dispersion. 
The points at which Ek- = are called Fermi points. 



The value of the Zeeman field at which new Fermi points 
appear represents a q uantum critical point for topological 
phase transition [l2l [l9| . For fi > 0, the system may 
support 0, 2 or 4 Fermi points along the fc z -axis as h 
is increased. For fx < which is the case illustrated in 
Fig. HI there can be either or 2 Fermi points [l3j]. For 
Xkp/Ep = 2, the critical Zeeman field is h c « 1.5Ep- 
For h < h c , the system is a topologically trivial gapped 
supcrfluid; for h > h c , the system possesses two Fermi 
points and represents a topologically nontrivial gapless 
superfluid. Note that the quasi-particle dispersion has 
been measured in recent experiments on ultracold Fermi 




— ■ , , , , , I i , , , , 

1 2 3 1 2 3 

k/kr k/kf 

FIG. 4: (color online). Quasi-particle dispersion spectrum 
iJk+ (solid lines) and -Ek- (dashed lines) shown in in units of 
Ef for k in the transverse plane (a, 6 = it/2) and along the 
2-axis (b, = 0) for \k F /E F = 2, h/E F = 2. 

Even though in our model the interaction has a contact 
s-wave form, due to the presence of the SO coupling, 
pairing can occur in both singlet and triplet channels 
[24| . The singlet pairing field between unlike spins can 
be calculated as: 



7 



(V'ktV'-ki) = -A < 



2/Cgkg) - 1 
4£ kQ 



1 



while the triplet pairing fields between like spins are given by: 

X(k y + ik x )(£k - h) 



(V'ktV'-kt) = A 
(V'kJ.V'-kJ.) = -A 



A(fc y - ifecXfk 



/o 



2/Cgk+) - 1 2/(£ k _) - 1 
2/(£ k _) - 1 



v /(^ 2 + A2^)^ + /i 2 A^ 



4£ k+ 

2/(ff k +) - 1 
4£ k+ 



4£ k _ 



(29) 

(30) 
(31) 



We plot in Fig.[5]the absolute value of the various pairing 
fields. One can see that the effect of the Zeeman field is 
to reduce the singlet pairing while enhancing the triplet 
pairing. 




0.25 



0.50 



0.5 1.5 2.5 

kjte 



FIG. 5: (color online). Pairing fields at unitarity (i.e., l/a 3 = 
0) for h = (top row), h = IEf (middle row), h = 2Ef 
(bottom row) and A/cf/Sf = 2. In each row, from left to 
right, we display |{V>k#-k4.)l> | (V'ktV'-kf) I and | (Vt^-H) | 
(all in units of Ef), respectively. 

From the single-particle Green function, one can imme- 
diately obtain the density of states which is an important 
quantity characterizing the nature of the quantum state. 
To this end, we invert Eq. (|23)l to get 



So (k, ioj m ) 



g(k,iu m ) 
/(+k, -iu) m ) 



/(k, iu m ) 

[g(-K -iuj m )Y 



(32) 

where g and / are both 2x2 matrices whose expressions 
can be found Appendix [3] The density of states are 
related to the g matrix as 





1 


— Im 




7T 


V 




1 


— Im 




7T 


V 



-^ 5ll (fc, w + ?:o+) 

k 

^322(fc,^ + i0+) 



(33) 
(34) 



where the expressions of gn and 1722 can be found in 
Eqs. (|A5[) and (| A7|) in the Appendix, respectively. 
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FIG. 6: (color online). Density of states p\ (a) and p_i (b) at 
unitarity and XIcf/Ef = 2. 

We show in Fig. [6] the density of states at different 
Zeeman fields. At h = 0, p\ and pi are identical, both 
exhibiting a large gapped region. As h increases, the 
gapped region shrinks. Furthermore, the density of states 
for the majority component (which is the spin down com- 
ponent for our choice of positive h) becomes more V- 
shaped near oj — 0. A V-shaped density of states is char- 
acteristic of many unconventional superconductors fea- 
turing non-s-wave pairing. Here it is due to the triplet 
pairing resulting from the SO coupling and enhanced by 
the Zeeman field. The density of states may be measured 
in experiment using the scheme proposed in Ref. 33j. 

Finally, let us consider spin structure factor [33, Hf| 
which is related to the dynamic spin susceptibility 
Xs (q, iy n )i which is the Fourier transformation of the 
spin-spin correlation function 



Xs (x,r) 



{T T 5n s (x,t) Sn s (0,0)} 



(35) 



where the spin density is given by Sns — — n\,. Us- 
ing the Nambu spinor notation, the spin density can be 
written as, 



Sns = -$ + 
6 2 



10 
0-100 
0-10 
1 



$ = -$+ [ 

1 I- 



I*. 



(36) 



8 



where Tj are the Pauli matrices describing the Nambu 
spinor degrees of freedom. Within the static Bogoliubov 

I 



approximation |36j, we have, 



Xs (q, iv n ) = j k BT ^ Tr {T* ® &zGo (k, iu! m ) f z ® <7 Z S (k + q, iuj m + iv n )} . 



(37) 



The zero-momentum dynamic spin structure factor is given by 



1 1 r i 

-f-lm [ XS (q = 0, 



1 — e —uj/k B T J]- 



o + )] , 



(38) 



r 



whose analytic expression is too complicated to show. 
However it can be explicitly shown that the dynamic spin 
structure factor vanishes when A = 0, i.e., in the absence 
of the SO coupling. Fig. [7fa) shows dynamic spin struc- 
ture factor at different Zeeman fields. At h = 0, Ss(0,u>) 
exhibits a broad peak. For finite h, an additional nar- 
rower peak appears at smaller energy. The corresponding 
static spin structure factor is given by 



S s (0) = / duiSsiO^), 

which is plotted in Fig. EJb) as a function of the SO cou- 
pling strength. The spin structure factor may be directly 
measured in experiment using the Bragg spectroscopic 
method 13711. 
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FIG. 7: (color online), (a) Zero temperature dynamic spin 
structure factor Ss(0, w) at unitarity and XkF / Ef = 2. (b) 
Static spin structure factor Ss (0) as functions of the SO cou- 
pling strength. 



V. OUTLOOK AND CONCLUSION 

In this paper, we provide a detailed theoretical descrip- 
tion of a two-component Fermi gas with Rashba SO cou- 
pling, under the framework of the functional path integral 
formalism. Many important and experimentally relevant 
physical quantities in both two- and many-body situa- 
tions — such as the effective mass of the two-body bound 
state, the chemical potential, gap parameter, pairing 



correlation, quasi-particle dispersion, density of states, 
spin structure factor, etc. — are calculated in a rather 
straightforward way using this method. The SO cou- 
pling in general favors the formation of Cooper pair and 
bound state. In particular, we have found that when a 
Rashba ring exists for the single-particle ground state, a 
two-body bound state exists regardless of the sign of the 
s-wave scattering length. 

For the sake of simplicity, we have only presented the 
zero-temperature results in this paper. The formalism 
and the relevant equations we have derived are neverthe- 
less valid for finite temperatures. In addition, we have 
focused on a uniform system. To take the trap into ac- 
count, one may adopt the local density approximation 
[H EH ; under which the trapped Fermi gas was treated 
as the sum of many small cells with a local chemical po- 
tential /i(r) = (j, — Vtrap(r). In our earlier work 9], we 
have shown that the important qualitative features of the 
system, like the mixed-spin pairing, anisotropy, enhanced 
pairing field, etc., are not affected by the presence of the 
trap. A self-consistent way to include the effects of the 
trapping potential is to use the Bogoliubov-de Gennes 
(BdG) formalism. For large number of atoms as used in 
most experiments, we expect that the local density ap- 
proximation should be quite accurate. However, for a 
system with relatively small number of atoms and tight 
trapping confinement, or in the case of vortex state where 
the order parameter varies significantly in a short length 
scale, the BdG approach will be more appropriate. 

One crucial feature of the system arising from the SO 
coupling is that the superfluid transition temperature is 
greatly enhanced [1(3] ■ The mean-field calculation which 
we have focused on in this work is not expected to pro- 
vide an accurate estimate of the transition temperature, 
particularly for system with strong interaction strength. 
As we have outlined in Sec. |H] of this paper, the Gaus- 
sian fluctuations on top of the mean-field level can be 
accounted for using the functional path integral formal- 
ism. This will lead to a much more accurate calculation 
of the transition temperature. We plan to address this 
issue in a future work. 
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Appendix A: Expressions of Single-Particle Green Function 

The single-particle Green function is given in Eq. f|32[) where the explicit expressions of the elements are 



with 



and 



M 



g(k, iu> m ) = [ioj m + £ k - ha z + \{k y a x - k x a y )} T>+/M, 
/(k,iw m ) - [-iA a y ]V_/M, 



(»W m - hf - £ 2 - A 2 - ( - Wm + h y _ ^2 _ A 2 _ ^2 ^2^2 ^2 _ ft 



V+ = 
V- -- 

In greater detail, we find that, 



where. 



(iw m + hf - t& - A 2 - \ 2 k\ A (k y + ik x ) [2 (£ k + ft)] 
A (k y - ik x ) [2 (& - ft)] (io; m - ft) 2 - ^ - Ag - A 2 fcj 



ftf-^-A 2 -^ 2 



-\{k v -ik x ) [2 + 
-A + ifc B ) [2 (a - ft)] (iu m + hf - £ 2 - A 2 , - X 2 k 2 ± 



5n(k,MJ m ) 3i2(k,iw m ) 
5*2 (k, -iw m ) 522 (k, zw m ) 



(Al) 

(A2) 
(A3) 

(A4) 



and 



M gil (k,iu) m ) = [ioj m + &-h] (iu) m + h) - £ - Ag - \ 2 k± +2\ 2 ki(&-h) 



Mg 12 (k, iuj m ) = A (fc 9 + ik x ) 



(iuj m + £ k ) 



A^ - \ z k ± 



Mg 22 (k,iuj m ) = [ioJ m +Zu + h} {ioj m - hf - £ - A 2 - X 2 k]_ + 2\ 2 k\ (& + h) , 



/(k, iu) m ) = ^ 



A(fcj,+jfe x )[2 (£ k -ft)] 
[iw m - hf - & - A 2 - \ 2 k\ 



(iuj m + ft) -t£ 



A 2 fc 2 



-A(fc y -i^)[2 (£ k + ft)] 



(A5) 
(A6) 
(A7) 

(A8) 
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